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I.  INTRODUCTION 


Me  have  recently  developed  a  one-dimensional ,  premixed,  laminar  flame 
code  that  considers  detailed  elementary  reactions  and  detailed  transport 
propertlesl.  As  a  test  of  this  code  we  elected  to  examine  the  ozone  flame, 
since  *ts  flame  speed  has  been  measured^  over  the  range  17  to  100  mole  % 


ozone.  The  reactions 
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As  expected,  the 

most  important 

reaction  for  our 

flame  studies 

proved  to  be  reaction 

(1).  Since  adiabatic  flame  temperatures  range 

from  1027  K  to  2677  K 

we  required  a 

rate  expression 

for  reaction  ( 

valid  over  this  temperature  range. 


In  196S  these  reactions  were  reviewed  in  detail  by  Johnston^.  He 
found  that  ki  (M  «  0j)  ■  9.94  x  10l4  exp( -22 . 72/RT)  cm-'  mole‘1  s-1  over 
the  temperature  range  200-1000  K .  In  1976  the  Leeds  group4  rccoranonded 
his  value.  (The  value  for  R  here  is  taken  as  1.9872  kcal  mole-1  K“l. 

To  convert  to  SI  units  note  that  4.184  joules  ■  1  calorie.) 


The  only  data  in  the  temperature  range  "v  600-1000  K  available  to 
Johnston  were  those  of  Jones  and  Davidson^.  AH  other  direct  measure¬ 
ments  of  reaction  (1)  used  by  Johnston  lay  in  the  range  303-5S9  K,  with 
most  of  these  taken  at  temperatures  less  than  400  K^. 
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cies”,  SSRDS-SBS20,  September  1968. 

4.  D.  L.  Baulah,  D.  D.  Drysdale,  J.  l^uxbury  and  S.  J.  Grant,  "Evaluated 
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1976,  Vol.  3. 
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Several  investigators,  all  using  shock  tubes,  found  indications  that 
the  Arrhenius  expression  for  reaction  (1)  given  by  Jones  and  Davidson8 
■sight  not  be  valid  at  higher  temperatures .  Over  a  combined  temperature 
range  of  1340-3300  K,  hrayb  and  Kiefer  and  Lutz'.8  found  their  experi¬ 
mental  results  were  consistent  with  values  of  k j /k2  that  were  much  lower 
than  those  predicted  using  the  Jones  and  Davidson  expression.  Michael^ 
used  Benson  and  Axworthy’slQ  value  for  kj  and  found  that  his  values  for 
kj  lay  progressively  lower  than  those  obtained  by  assuming  the  Jones  and 
Davidson  expression  to  be  valid  over  the  range  "v  1000-1400  K.  Finally 
Center  and  KungH  showed  experimental  evidence  that  this  rate  coefficient 
is  about  a  factor  of  two  smaller  than  the  value  obtained  by  assuming  that 
Johnston's  expression  is  valid  at  3000  k. 

The  expressions  of  Johnston  and  of  Jones  and  Davidson  describe  the 
temperature  dependence  of  ki  over  a  limited  range.  Sufficient  informa¬ 
tion  seemed  to  be  available  to  arrive  at  a  better  description  of  the 
temperature  dependence  of  reaction  (1).  The  data  of  Michael^  and  the 
high  temperature  data  of  Center  and  Ming*-  proved  to  be  in  a  form  that 
made  this  possible.  \  complete  listing  of  the  low  temperature  data  com¬ 
piled  by  Johnson8,  essentially  his  Table  XVII 1,  and  that  of  Michael^  is 
listed  in  Appendix  A.  The  data  ot  Center  and  kungl-  are  listed  in 
Appendix  B. 


II.  ANALYSIS 

A.  General 

he  desired  values  for  kj  extending  from  300-3000  k  so  that  a  valid 
fit  (i.e.  Arrhenius  description)  could  be  obtained.  We  should  like  to 
use  as  much  of  Michael's  data  and  of  Center  and  Ming's  data  as  possible 
to  extend  the  temperature  range  covered  in  Johnston's  review.  To  this 
end  we  first  examined  the  dependence  of  Michael's  data  for  k(  upon  the 
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FA,  1967,  pp  67-76. 
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value  of  k2  and  then  determined  appropriate  high  temperature  values  for 
ki  from  the  data  of  Center  and  Kung. 

In  evaluating  both  the  data  of  Michael  and  of  Center  and  Kung  we 
have  followed  Johnston  and  have  converted  each  of  their  measurements  to 
equivalent  ozone.  To  do  so  we  have  used  kj  (M»Ar)/ki(M«Kr)  «  1.2S 
(ref.  9)  and  kj  (M-OjJ/kj (M«Ar)  ■  4.0  (ref.  3)  and  have  assumed  that 
these  ratios  are  independent  of  temperature.  Should  quantitative  third 
body  efficiencies  as  a  function  of  temperature  become  available  and  be 
found  to  differ  significantly  from  the  values  we  have  employed,  the 
analysis  presented  below  would  have  to  be  repeated  using  appropriately 
adjusted  data. 


B.  Michael 's  Data 
q 

Michael  operated  his  shock  tube  at  low  initial  pressures  0.04  - 
0.2  atmospheres)  and  large  amounts  (^  9S\)  of  krypton  as  a  diluent  and 
so  the  only  reactions  of  importance  are  (1)  and  (2).  This  leads  to  the 
following  rate  equations: 


d(03)/dt  -  -  kjlOjUM)  -  k2[03][0]  (4) 

and 


dJOj/dt  -  ♦  k j lOj] I M)  -  k2l03]|0).  (S) 


When  the  oxygen  atom  concentration  reaches  the  steady  state;  i.e., 
d[0]/dt  «  0,  we  can  write 


d(03]/dt  -  -  2  kj(03](M). 


(6) 


Reactions  (1)  and  (2)  were  modeled  by  Michael  and  the  logarithm  of  the 
calculated  ozone  concentration  was  plotted  against  time.  He  found  no 
significant  deviation  from  first  order  kinetics,  and  kl  was  described  by 
equation  (6).  However,  he  had  used  Benson  and  Axworthy's  value*0  for 
k2  •  3.37  x  10*5  cxp(-S ,70/RT)  cm3  mole**  s-*.  Cater  expressions  of 
Johnston3,  k^  -  1.20  x  10*3  exp(-4 . 79/RT)  cm3  mole-*  s**  (200-1000  K) , 
and  the  Leeds  group4,  k2  ■  S.2  x  10*2  exp( -4 . 1S/RT)  cm3  mole-1  s**  (200- 
S00  K)  lead  to  sma 1 1 cr  values  for  k2.  This  means  that  the  time  for  the 
oxygen  atom  concentration  to  reach  steady  state  will  take  longer  than  in 
Michael's  calculations.  The  question  is:  how  much  longer? 


We  have  numerically  intergrated  equations  (4)  and  (5)  for  Michael's 
conditions  using  the  smallest  value4  for  k2  and  found  that  d[0)/dt  3  0 
is  still  valid  in  the  time  frame  of  Michael’s  experiment,  typically  50 
to  350  usee.  (Rate  coefficients  derived  using  equation  (6)  will  tend  to 
be  systematically  low  but  the  error  is  less  than  1 0\) .  Note  that  we 
must  assume  that  the  expression  for  k2  is  valid  up  to  s.  1400  K.  This 
approach  appears  to  be  reasonable  since  Michael's  data  are  not  very  sen¬ 
sitive  to  the  precise  value  of  kj. 

In  summary  we  find  that  the  later,  lower  values  of  k2  do  not  signi¬ 
ficantly  affect  Michael's  results. 

C .  Center  ji nd  kung '  s  Da t.i 

Center  and  kungll  operated  their  shock  tube  at  very  low  initial 
pressures  ps  7  x  10"*  to  3  x  10‘-  atmospheres)  with  a  high  ( 95 -99\) 
argon  diluent.  The  time  scale  of  these  measurements  is  tens  of  micro¬ 
seconds,  much  shorter  than  the  time  scale  in  Michael's  experiment.  A 
simple  analytical  unfolding  of  their  ozone  concentration  relaxation  time 
observations  could  be  made  provided  that  |0)  0  or  more  precisely  that 

k2  10)  <<  ki  |M).  Under  these  conditions 


d[0.j/dt  3  -kj  [M]  (03J.  (7) 


Numerical  integration  of  equations  (4)  and  (5)  for  their  experi¬ 
mental  conditions  shows  that  as  one  progresses  from  2000  k  to  lower  temp¬ 
eratures  the  more  important  the  term  involving  atomic  oxygen  becomes 
(see  liquation  4),  and  that  for  temperatures  greater  than  2000  k.  the  loss 
rate  of  ozone  can  be  approximately  described  by  equation  (7).  In  this 
manner  fourteen  points  of  Center  and  kung*-  were  determined  to  be  valid 
representations  of  ki  and  arc  1 i sted  in  Table  I .  A  complete  listing  of 
their  data  and  a  more  detailed  analysis  of  their  data  may  be  found  in 
Appendix  B. 

The  question  remains:  what  is  the  sensitivity  of  these  values  of 
kj  as  listed  in  Table  I  to  the  value  of  k2?  It  can  be  seen  in  equation 
(4)  that  the  larger  the  value  of  k2  the  more  important  that  term  becomes. 
Here  we  have  used  Mampson's*'  expression  for  kj  *  1.14  x  1013  cxp(-4.57/ 
RT)  cm3  mole* l  s‘*  (200-1000  k) .  Again  this  expression  for  k2  has  been 
assumed  to  be  valid  at  the  higher  temperatures.  To  test  the  sensitivity 
of  the  values  of  kj  on  the  value  of  k’  we  have  arbitrarily  multiplied 
liampson's  expression  for  k2  by  two  for  the  2041  K  case  (see  Table  I). 

We  have  found  that  the  computed  ozone  relaxation  time  (i.c.,  kj)  changes 
by  about  10\.  For  the  higher  temperatures  the  change  is  less. 

13.  F.  F.  Hjnpean,  ED.t  •/.  rhua.  and  ('Turn.  Ref.  Data  £,  267-308  (1973). 


TABLE  I.  HIGH  TEMPERATURE  VALUES  OF  kj  FROM  CENTER  AND  KUNG 


T(K) 

k;  (M  ■  Ar) 

2041 

5.31(11)* 

2128 

9.48(11) 

2273 

1.54(12) 

2353 

1.76(12) 

2439 

1.18(12) 

2500 

9.98(11) 

2564 

1.08(12) 

2564 

1.30(12) 

2632 

1.37(12) 

2667 

1.62(12) 

2703 

1.84(12) 

2778 

2.25(12) 

2857 

2.24(12) 

2941 

1.90(12) 

Read  ue  5.31  x  10“  am'  mole' "  o' “ . 


D.  The  Least  Squares  Fit 


The  data  can  be  divided  into  four  groups  by  temperature:  the  lower 
temperature  values-*  (303-359  KJ ;  the  Jones  and  Davidson^***  values,  (769- 
910  K);  Michael's9  values  (971-1384  K) ;  and  Center  and  Kung'sH  values, 
(2041-2941  K) .  So  that  the  preponderance  of  low  temperature  data  points 
would  not  bias  the  fit,  each  group  of  data  was  weighted  equally.  For 
example,  we  used  14  values  of  kj  from  Center  and  Kung  (see  Table  I)  and 
in  this  group  each  value  was  weighted  by  1/14.  To  avoid  biasing  the  fit 
toward  the  higher  absolute  values  of  k j ,  errors  between  the  least  squares 
fit  and  the  data  points  were  all  measured  in  a  relative  sense.  This  is 
accomplished  by  using  a  weight  for  each  value  of  kj  equal  to  the  square 
of  the  reciprocal  of  that  value.  This  weight  is  in  addition  to  the 
weighting  discussed  above.  The  logarithmic  method  with  proper  transfor¬ 
mation  of  the  weights*-*  was  used  to  fit  the  data  points  to  the  two  para¬ 
meter  Arrhenius  form:  A  exp(-E/RT),  where  R=1  .9872  x  10"  -*  kcal  mole'*  K~ * 
and  T  is  the  temperature.  The  value  k(  (M=03 )  *  4.31  x  10*-*  exp(-22.2/ 

RT)  cm-'  mole*  s**  was  found.  The  number  of  digits  carried  is  a  measure 
of  the  precision  of  the  fit  and  not  of  the  accuracy  of  the  experiments. 

A  listing  of  the  least  squares  fitting  routine  coded  in  FORTRAN  for  use 
on  a  CDC - 7&00  can  be  found  in  Appendix  C. 

Table  II  shows  the  average  error  for  the  points  in  each  of  the  four 
data  groups  relative  to  our  recommended  fit.  The  average  error  for  all 
the  data  points  is  also  shown.  The  fact  that  the  average  relative  error 
for  each  data  group  is  comparable  to  the  error  in  the  overall  fit  is 
another  indication  that  the  fit  is  not  dominated  by  any  single  data  group. 
The  site  of  the  errors  can  be  traced  to  the  large  amount  of  scatter  in 
the  data. 

To  check  the  consistency  of  these  independent  data  and  the  robust¬ 
ness  of  the  derived  two  parameter  Arrhenius  expression,  fits  were  also 
made  by  deleting  one  group  of  data  at  a  time.  These  results  are  given 
in  Table  III.  The  evaluation  of  these  two  parameters  is  seen  to  be  in¬ 
dependent  of  any  single  data  group.  Specifically,  if  the  highest  tem¬ 
perature  data  of  Center  and  Kung  are  deleted,  the  subsequent  expression 
for  kj  extrapolated  to  2500  K  yields  values  only  23  greater  than  all  the 
data  combined .  Similarly  if  Michael's  data  are  deleted  the  subsequent 
interpolated  value  of  this  expression  for  kj  is  173  greater  than  all  the 
data  combined. 

Separate  fits  were  made  to  each  of  the  four  data  groups.  In  some 
cases  the  parameters  so  obtained  differed  markedly  from  those  obtained 
using  all  the  data.  This  comes  as  no  surprise  since  each  of  the  four 
groups  encompass  a  rather  limited  temperature  range,  while  taken  together 
they  span  a  temperature  range  of  about  a  factor  of  ten. 


14.  F.  J.  Cvetanovic  and  D.  I.  Singleton,  Intemat.  J.  of  Chem.  Kin.  9, 
481-488,  1977.  See  also,  F.  J.  Cvetanovic  and  D.  L.  Singleton, 
Intemat.  J.  of  Chem.  Kin  9,  1007-1009 ,  1977. 
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TABLE  II.  COMPARISON  OF  DATA  FROM  EACH  GROUP 
WITH  OUR  RECOMMENDED  FIT 


Data  Group 

Average  Relative  Error 

Low  Teap 

251 

Jones  6  Davidson 

43\ 

Michael 

41% 

Center  6  Kung 

25% 

All  Data 

33% 

TABLE  III.  LEAST  SQUARES  FIT  FOR 
EXCLUDING  DATA  GROUPS  ONE  BY 

ALL 

ONE 

DATA 

All  Data  Except 

Log  A 

E (kcal/mole) 

Low  Trap 

14.39 

20.6 

Jones  6  Davidson 

14.60 

22.2 

Michael 

14.70 

22.3 

Center  6  Kung 

14.64 

22.2 

All  Data 

14.63 

22.2 

11 


A  three  parameter  fit  of  the  data  was  found  to  yield  the  expression 
kj  *  5.31  x  1016  t-0.61  exp(-23. 1/RT)  cm3  mole*!  s-1.  The  activation 
energy  is  about  the  same  as  the  two  parameter  case  and  the  average  error 
for  all  the  data  in  this  case  is  31\,  slightly  smaller  than  the  two  para¬ 
meter  case.  The  three  parameter  fit  is  sensitive  to  data  deletion  and 
thus  we  feel  that  the  data  are  not  sufficiently  precise  to  support  such 
a  fit.  The  two  parameter  fit  is  essentially  just  as  accurate. 

Figure  la,  b,  and  c  show  a  plot  of  all  the  data  together  with  both 
the  two  parameter  fit  (solid  line)  and  the  three  parameter  fit  (dashed 
line)  to  these  data.  Since  the  vertical  scale  spans  about  15  powers  of 
ten,  the  plot  has  been  segmented  into  three  sections. 

For  T  3000  K  the  two  parameter  expression  we  have  derived  leads 
to  values  of  ki  that  arc  about  a  factor  of  two  lower  than  those  obtained 
by  assuming  that  the  Johnston  expression  is  valid.  Because  the  Jones 
and  Davidson  data  constituted  the  highest  temperature  available  to  him, 
Johnston  necessarily  weighted  his  fit  heavily  toward  these  data.  However, 
as  can  be  seen  in  the  figure,  our  expression  yields  values  for  kj  that 
lie  below  al 1  the  Jones  and  Davidson  data.  We  do  not  know  if  the  Jones 
and  Davidson  measurements  are  systematically  too  high^,  0r  if  there  is 
actually  some  pre-exponential  temperature  dependence  in  the  function 
representing  kj  that  cannot  yet  be  discerned  because  the  data  arc  not 
sufficiently  accurate. 
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For  a  diecueeion  of  some  data  reduction  problems,  see  H.  B.  Palmer 
I'onbusticm  and  Flane  11,  120-124,  1967. 


Figure  la. 
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*  HICHREL  \ 
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1000/T(K) 

High  Temperature  Data  for  the  Uniroolecular  Ozone  Decomposition 
Reaction.  The  data  of  Center  and  Kung  arc  taken  from  Table  1, 
the  data  of  Michael  from  reference  9  of  text  and  the  data  of 
Jones  and  Davidson  from  reference  3  of  text.  The  data  in  the 
figure  have  been  adjusted  for  M  *  O3.  The  solid  and  dashed 
lines  show  respectively  the  best  two  and  three  parameter  fits 
to  all  the  data. 
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LOG  KKCM3/MOLE-SEC) 


1.3  IX  1.7  1.9  8.1  2.3  8.6 


1000/T( K  \ 


ligurc  lb.  Intermediate  rcaperature  Range  Data  for  the  Unimolccular  Ozone 
Decomposition  Reaction.  These  data  are  taken  from  reference 
3  of  text .  Table  18 . 


LOG  K 1 l CM3/M0LE - SEC ) 


1000/T(K) 


Figure  1c.  l.ow  Temperature  Range  Pata  for  the  Unimolccular  O:one  Decom¬ 
position  Reaction.  These  data  are  taken  fro*  reference  3  of 
text  ,  Table  18. 
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APPENDIX  A.  DATA  EXCLUSIVE  OF  CENTER  6  RUNG 


The  purpose  for  this  appendix  is  to  list  ail  the  data  used  for  this 
evaluation  of  kj  so  that,  together  with  Table  I,  it  would  be  available 
in  one  location  for  future  use. 

The  bulk  of  Table  1-A  that  comprises  this  appendix  is  taken  directly 
from  Johnston  (reference  3  of  text) .  The  only  change  from  his  Table  18 
is  the  deletion  of  the  six  data  points  for  which  Axworthy  added  H2O2  or 
H2O  to  the  gaseous  mixture.  The  references  for  this  table  are  to  be 
found  in  reference  3  (of  text).  Following  Johnston,  these  data  are 
listed  for  equivalent  ozone  (i.e.,  M  •  O3),  and  we  have  converted  these 
to  cm^/moles  using  the  conversion  6.023  x  10^3  particles/mole. 

Michael's  data  are  taken  from  reference  9  (in  text)  where  equivalent 
ozone  has  been  found  using 


kj  (M  -  ArJ/kj  (M  »  Kr)  -  1.25 


and 


kj  (M  -  03)/kj  (M  -  Ar)  -  4.0. 


Tabic  1-A  is  comprised  of  three  coliaans.  The  first  lists  degrees 
Kelvin,  the  second  column  lists  the  reciprocal  of  the  first  multiplied 
by  1000.  The  third  column  lists  the  logarithm  to  the  base  ten  of  the 
rate  coefficient  for  reaction  (1)  in  units  of  cmVmolc-sec •  All  these 
values  are  given  for  ozone  as  the  third  body.  To  convert  to  other  third 
bodies  see  references  3  and  9  (of  text). 
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APPENDIX  B.  DATA  OF  CENTER  AND  RUNG 


The  purpose  of  this  appendix  is  two  fold.  First  and  foremost  to 
serve  as  a  repository*  for  the  actual  tabular  data  upon  which  Center  and 
Rung  have  published  their  paper  (see  Figure  4  of  reference  11  in  text). 
The  second  purpose  lies  in  providing  more  detail  to  the  analysis  given 
this  data. 

Center  and  Kung's  data  (reference  12)  are  given  in  Table  B-I.  The 
terms  employed  are  defined  in  the  Glossary  for  this  appendix.  In  ana¬ 
lyzing  this  data  set,  the  main  problem  lies  in  finding  a  value  for  kj 
given  the  measured  relaxation  time,  Tp. 

Provided  k2(0]  <<  k j [M]  in  equation  4  (in  text)  we  can  write 
d(0j]/dt  *  -kj (MJ (Oj] .  Since  their  measurements  are  made  after  the  shock 
we  have 


l"l  '  -2  I  Wr  •  4''o5/"a5]  • 


(B-l) 


where  following  Johnston  (reference  3)  we  have  taken  kj  (M  «  O3)  «  4k  j 
(M  •  Ar) .  Since  9S\  or  more  of  the  gas  is  argon  it  is  a  good  approxi¬ 
mation  to  take  [Ml  ■  constant.  Then  equation  7  (in  text)  can  be  readily 
integrated,  to  wit 


103]  "  l O3 ) 0  «*P  t-kjlMjt). 


(B-2) 


The  experiment  consists  of  optically  monitoring  the  ozone  concentration 
and  noting  the  time,  Tp,  that  the  post  shock  intensity  equals  the  pre¬ 
shock  intensity.  This  corresponds  to  a  change  in  the  O3  concentration 
equal  to  02/01 .  Then 


fn(p j/p^)  •  -k j I M 1  t  . 


(B-3) 


We  want  an  expression  for  ki  in  terms  of  experimentally  measurable  para 
meters.  Using  equation  (B-l),  equation  (B-3)  can  be  rewritten  as 


*Dr.  Kung  has  graoiously  consented  that  this  be  done  since  these  data 
do  not  appear  anjvhere  else. 


2  3 


£n  Ip,/pj) 


(B-4) 


(»,/»,)  (»,.p)  (VAr/«Ar  .  .Y0j/»0j) 


An  expression  for  Pj/Pj  can  be  found  from  the  ideal  gas  law: 


Vi  ■  lRTi  <  W 


(B-S) 


and  measuring  r  in  usee  we  find 


k,  (M  -  Ar)  • 


10  <n  (o,/p,)  «t,  (rAr/«Ar  .  >„  /«„  ) 

(Pi’p1  (W  •  4  V/V 


cm'  mole  *s  (B-6) 


Before  accepting  any  value  of  kj  derived  from  equation  B-b,  we  must 
verify  the  basic  assumption  that  k’  |0]  <<  kj  |MJ.  To  do  this,  equations 
(4)  and  IS)  (in  text)  are  integrated  numerically.  Because  of  the  uncer¬ 
tainty  in  using  expressions  for  the  value  of  k^  at  the  higher  tempera¬ 
tures.  a  value  twice  as  large  as  the  largest  recommended  value  (reference 
13  of  text)  has  been  used  in  the  numerical  integration.  Those  integra¬ 
tions  that  gave  relaxation  times  within  10\  of  the  measured  relaxation 
times  were  accepted  as  meeting  the  criterion  kj  |M]  <<  k 2  1 0 J . 

For  example  at  Tj  »  2041  K,  the  measured  relaxation  time  is  13.4  usee. 
If  Kampson’s  value  (reference  13  of  text)  for  k^  is  used,  wc  calculated 
a  12.7  usee  relaxation  time.  If  twice  Hampson’s  value  is  used,  the  re¬ 
laxation  time  is  found  to  be  12.2  usee.  Since  the  term  ki  (0J  in  equa¬ 
tion  4  (of  text)  becomes  the  more  important  the  lower  the  temperature, 
the  datum  at  Tj  ■  2041  K  is  the  lowest  temperature  we  have  accepted. 
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TABLE  B- 1 .  THE  DATA  COLLECTED  BY  CENTER  AND  KUNG 


Run 

V 

a 

s 

CMI/USJ 

Pi 

(torrj 

p2/cl 

Vi  ’jpi 

(•taos-uscc I  (ataos-wsec) 

1000/T 

IK'1) 

686 

6 

1.27 

0.74 

5.645 

4.5  (-2)- 

5.1  (-2) 

0.65 

68' 

6 

1.55 

0.758 

5.86 

7.7  (-2) 

5.6  (-5) 

0.47 

688 

5 

1.4S 

0.725 

5.8 

1.2  (-5)** 

9.0  1-5) 

0.51S 

690 

S 

1.515 

0.70 

5.69 

5 .8S( *2) 

2.5  (-2) 

0.60 

691 

5 

1.50 

0.682 

5.67 

5.9  (-2) 

2.5  (-2) 

0.61 

69.’ 

5 

1.28 

0.675 

5.65 

4.2  (-2) 

2.8  (-2) 

0.62 

694 

S 

1 .27 

0.6SS 

5.645 

4.8  (-2) 

5.0S( -2) 

0.65 

69S 

s 

1.725 

0.650 

5.97 

4.6  (-5) 

1.9  (-5) 

0.575 

698 

5 

1.42 

0.620 

5.77 

2.0  (-2) 

0.55 

699 

S 

1.21 

0.610 

5.58 

0.70 

S.S  (-2) 

0.67 

700 

s 

1.12 

0.605 

5.47 

0.145 

0.155 

0.74 

701 

s 

1.06 

1.68 

5.58 

0.255 

0.205 

0.78 

70  2 

& 

1 .04 

2.05 

5.55 

0.27 

0.25 

0.80 

705 

s 

1.095 

2.41 

5.45 

0.247 

0.145 

0.75 

710 

1  .59 

0.775 

5.75 

8.2  (-5) 

S.S  (-5) 

0.40 

7 1 S 

y 

0.96 

2.72 

5. IS 

0.91 

0.74 

0.92 

716 

y 

1  .01 

2.57 

5.21 

0.61 

0.42 

0.86 

'44 

& 

1.57 

0.767 

5.75 

2. 2S( -2) 

1 .  S2 ( -2) 

0.565 

74S 

& 

1.65 

0.7S5 

5.95 

6.5  1-5) 

2.851-5) 

0.41 

'46 

s 

1.69 

0.425 

5.95 

6.9  (-5) 

2.5  (-5) 

0.59 

747 

s 

1.27 

0.455 

5.65 

4.5  1-2) 

5.2  (-2) 

0.65 

748 

y 

1.42 

0.88 

5.65 

1 . 55( -2) 

1.0  (-2) 

0.49 

750 

y 

1.615 

0.86 

5.77 

6.5  (-5) 

5.1  (-5) 

0.59 

7  S 1 

y 

1.14 

1 .68 

5.59 

0. 14 

0.11 

0.72 

75: 

2 

1.175 

1  .65 

5.45 

9.6  (-2) 

7.2  (-2) 

0.70 

'55 

j 

* 

1.27 

1  .60 

5.55 

S.l  (-2) 

5.1  (-2) 

0.60 

7S4 

1.5 

1.255 

5.05 

5.4 

S.S  (-2) 

5.S  (-2) 

0.61 

7S7 

1.5 

0.99S 

2.82 

5.1 

0.62 

0.57 

0.88 

758 

1 .2 

0.955 

2.75 

5.1 

0.65 

0.S7 

0.88 

7S9 

1.2 

1.205 

1.58 

5.55 

7.7  (-2) 

S.S  (-2) 

0.655 

760 

1.2 

1.255 

2.62 

5.59 

4.2  (-2) 

4.2  (-2) 

0.625 

7&: 

1.2 

1.05 

2.51 

5.15 

0.67 

0.40 

0.855 

765 

1.2 

1.12 

2.48 

5.25 

0.246 

0.158 

0.755 

'66 

1.0 

1.5 

2.42 

5.44 

5.2  (-2) 

2.4  (-2) 

0.575 

767 

4.S 

1.815 

0.4 

4.01 

5.9S(-5) 

1.42(-5) 

0.54 

768 

4.S 

1.615 

0.582 

5.91 

4.5  (-5) 

5.4  (-5) 

0.425 

'69 

4 . 5 

1 .60 

0.578 

5.90 

4.9  (-5) 

5.0  (-5) 

0.44 

771 

4.5 

1.715 

0 . 580 

5.97 

S.S  (-5) 

2. 1 2 ( - 5 ) 

0.58 

774 

4 .  S 

1.77 

0.55 

5.99 

5. 5S( -5) 

1.7  (-5) 

0.56 

77S 

4.5 

1.745 

0.557 

5.98 

4.1  (-5) 

1.87(-5) 

0.37 

776 

4.5 

1 .79 

0.54 

4.00 

5 . 55( -5) 

1 . 56( -5) 

0.35 

'b: 

S 

1.41 

0.65 

5. '8 

1.6  (-2) 

1.2  (-2) 

0.54 

785 

5 

1.425 

0.615 

5.78 

1 . 191-2) 

1.1  (-2) 

0.525 

'84 

5/4 

0.99 

5.02 

5.09 

0.75 

.66 

0.89 

785 

5/4 

1.0 

2.96 

5.1 

0.75 

.57 

0.88 

'86 

5/4 

1.0 

2.84 

5.1 

0.69 

.57 

0.88 

'87 

5/4 

1  .06 

2.75 

5.18 

0.5SS 

.27$ 

0.80 

'89 

5/4 

1.05 

2.60 

5.17 

0.45 

.51 

0.81 

'90 

5/4 

0.97 

2.55 

5.06 

0.97 

.85 

0.91 

'91 

5/4 

1.055 

2.47 

5.15 

0.45 

.37 

0.85 

'88 

5/4 

1 .05 

2.68 

5.15 

0.46 

.4 

0.835 

4.&1-2)  it  4.$  x  id 
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GLOSSARY  FOR  APPLNDIX  B 


Xy  j  Initial  mole  fraction  of  0^. 
u^  Shock  speed  (mm  (1sec  *)■ 

Pj  Initial  pressure  (torr). 

.j  Initial  density  (gm/cm'). 

*  ,  Po>t -shock  density  (gm/cm'). 


t  Relaxation  time  measured  by  Center  and  Kung  (asec).  This  is  the 
time  in  which  the  0.  concentration  changes  by  a  factor  of  »  ,/c j . 

Relaxation  time  predicted  by  Johnston  recommended  rates  (^sec). 


Pj  Relaxation  time  multiplied  by  initial  pressure  (,<sec  -  atm). 

Tj  Initial  temperature  (k). 

Y^r  Mass  fraction  of  argon. 

Y  Mass  fraction  of  ozone. 

X 

l»  Molecular  weight  of  argon  (40  gm  mole). 

\  T 


k 


APPENDIX  C.  LEAST  SQUARES  FITTING  PROGRAM 

Below  is  a  listing  of  the  subroutine  used  to  compute  a  least  square 
fit  to  ki  using  the  data  listed  in  Appendix  A  and  Table  I.  The  program 
can  be  easily  modified  to  produce  three  parameter  fits  of  the  form 
aT*>  exp  (C/T). 
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